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Abstract 

A nonperturbative method for the solution of quantum field theories is described in the context 
of quantum electrodynamics and applied to the calculation of the electron's anomalous magnetic 
moment. The method is based on light-front quantization and Pauli-Villars regularization. The 
application to light-front QED is intended as a test of the methods in a gauge theory, as a precursor 
to possible methods for the nonperturbative solution of quantum chromodynamics. The electron 
state is truncated to include at most two photons and no positrons in the Fock basis, and the wave 
functions of the dressed state are used to compute the electrons's anomalous magnetic moment. A 
choice of regularization that preserves the chiral symmetry of the massless limit is critical for the 
success of the calculation. 

PACS numbers: 12.38.Lg, 11.15.Tk, ll.lO.Gh, ll.lO.Ef 



Contributed to the proceedings of QCD@Work 2010, the international workshop on QCD Theory and 
Experiment, Martina Franca, Italy, June 20-23, 2010. 



1 



I. INTRODUCTION 



The purpose of this work is to explore a nonperturbative method that can be used to solve 
for the bound states of quantum field theories, in particular QCD. The problem is notoriously 
difficult, and there are only a few approaches. These include lattice gauge theory the 
transverse lattice Q, Dyson-Schwinger equations [3], Bethe-Salpeter equation, similarity 
transformations combined with construction of effective fields jjl, light-front Hamiltonians 
with either standard [5i] or sector-dependent parameterizations |6|-l8|. 

We use the light-front Hamiltonian approach with Pauli-Villars (PV) jof regularizaton 
and standard parameterization, where the bare parameters of the Lagrangian do not depend 
on the Fock sector. This means that we use Fock states - the states with definite particle 
number and definite momentum for each particle - as the basis for the expansion of eigen- 
states. The coefficients in such an expansion are the wave functions for each possible set 
of constituent particles. These functions describe the distribution of internal momentum 
among the constituents. Such an expansion is infinite, and we truncate the expansion to 
have a calculation of finite size. 

The wave functions are determined by a coupled set of integral equations which are 
obtained from the bound-state eigenvalue problem of the theory. Each bound state is an 
eigenstate of the field-theoretic Hamiltonian, and projections of this eigenproblem onto in- 
dividual Fock states yields these coupled equations. Each equation is a relativistic analog 
of the momentum-space Schrodinger equation, but with terms that couple the equation to 
other wave functions that represent different sets of constituents, perhaps one gluon more 
or less or a quark-antiquark pair in place of a gluon or vice-versa. 

The solution of such equations, in general, requires numerical techniques. The equations 
are converted to a matrix eigenvalue problem by some discretization of the integrals or by 
a function expansion for the wave functions. The matrix is usually large and not diagonal- 
izable by standard techniques; instead, one or some of the eigenvalues and eigenvectors are 
extracted by the iterative Lanczos process. The eigenvector of the matrix yields the wave 
functions, and from these can be calculated the properties of the eigenstate, by considering 
expectation values of physical observables. 



We work with light-cone coordinates [10|, lU 



chosen in order to have well-defined Fock- 
state expansions and a simple vacuum. The time coordinate is x"*" = t + z and the space 
coordinates are x_ = [x~,x±), with x~ = t — z and x± = {x,y). The light-cone energy is 
p~ = E — pz, and the three-momentum is p = (p^, p±), with p"*" = E + pz and p± = [px, Py)- 



The mass-shell condition p^ 



m 



becomes p 



The simple vacuum follows from 
^ = ^y m'^ + pI + p^ +Pz > 0. 



the positivity of the plus component of the momentum: p 

To regulate a theory, we use the Pauli-Villars technique [9| . The basic idea is to subtract 
from each integral a contribution of the same form but of a PV particle with a much larger 
mass. This can be done by adding negative metric particles to the Lagrangian. For example, 
for free scalars a Lagrangian of the form 



1 1 



1 



2 i2 



generates a contribution from an internal line of a Feynman diagram in the form 

1 1 



p"^ - ^il p'^-fili 



d p, 



[1.2] 
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FIG. 1. The Fock-state expansion of the dressed-electron eigenstate. 

which has the necessary subtraction. A particular advantage of PV regularization is preser- 
vation of at least some symmetries; in particular, it is automatically relativistically covariant. 

II. APPLICATION TO QED 

The method is not mature enough to apply to QCD, so as a test in a gauge theory, 
we consider light-front QED and specifically the eigenstate of the dressed electron and 



its anomalous moment |12Hl5|. From the PV regulated light-front QED Lagrangian, we 



construct the Hamiltonian V and solve the mass eigenvalue problem V \P) = ^\P) in 
the approximation that the electron eigenstate is a truncated Fock-state expansion with at 



most two photons and no positrons ^5|. From this approximate eigenstate, we compute 
the anomalous magnetic moment from the spin-flip matrix element of the electromagnetic 
current J+ (l6j . 

Schematically, the electron Fock-state expansion can be written 

I electron) = / ^e|e) + / ?/'e^|e7) + / Ve77|e77) + / V'eee+|ee+) H (2.1) 



This is represented graphically in Fig. [T] It satisfies the eigenvalue problem 

i^LC [electron) = {K + Vqed) |electron) = M^jelectron, (2.2) 
where Vqed is the potential-energy operator, represented graphically in Fig. [2l The projec- 



Vqed = • + • *=^^— ' + 





FIG. 2. The potential-energy operator of light-front QED in a graphical representation. 

tions of this eigenvalue problem onto each individual Fock state produces coupled equations 
for the Fock-state wave functions, which are schematically represented in Fig. |3l The first 
graphical equation in Fig. |3] is a projection onto the one-electron Fock state, which can be 
written as 

mllPe + j dk^Ve^^e{k^)i^e^{k^) = M^e- (2.3) 
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FIG. 3. A graphical representation of the coupled equations for the Fock-state wave functions of 
the dressed electron eigenstate. 



It includes absorption of the photon from the one-electron/one-photon state. The second 
equation can be expressed as 



E 



^+^^^' ^e7 + J dk^Ve^^e{k^)^e-yy{k^) = MVe7- (2-4) 



This includes photon emission by the bare electron and photon absorption from the one- 
electron/two-photon state. The third equation is the analogous one for the three-body 
sector. 

The first and third equations of the coupled system can be solved for the bare-electron 
amplitudes and one-electron/two-photon wave functions, respectively, in terms of the one- 
electron/one-photon wave functions. Substitution of these solutions into the second integral 
equation yields a reduced integral eigenvalue problem in the one-electron/one-photon sector. 
A diagrammatic representation is given in Fig. |H 

Solution of the resulting integral equations yields a as a function of mo and the PV 
masses. Then for given values of PV masses, we can seek the value of mo for which a takes 
the standard physical value e'^/Air. The equations must first be solved for M = 0, with the 



coupling strength parameters adjusted to yield mo = [15[. The solution of the integral 
equations requires numerical techniques [l5[. The integrals are discretized via quadrature 
rules, and the equations are thereby converted to a matrix eigenvalue problem, which is 
solved by iteration. 
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III. RESULTS 

From the solutions to the eigenvalue problems, we compute the anomalous moment at 
fixed PV masses and fixed numerical resolution. We then study the behavior first as a 
function of the numerical resolution, which requires extrapolation, and then as a function 
of PV masses. The numerical resolution is marked by two parameters, K and N±, which 
control the number of quadrature points used in the longitudinal and transverse directions. 



The numerical convergence and extrapolation are illustrated in [15 . 
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FIG. 4. Diagrammatic representation of effective integral equation in the one-electron/one-photon 
sector. The first term represents the kinetic energy; the second and third, the two time orderings 
of photon absorption and emission; and the fourth, the self-energy contribution. 



The results of the extrapolations are plotted in Fig. [H Each value is close to the stan- 
dard Schwinger result of a/27r and independent of /ii, to within numerical error. The 
results with only the two-photon self-energy contribution are actually better than the full 
two-photon results. This discrepancy should be due to the absence of electron-positron con- 
tributions, which are of the same order in a as the two-photon contributions; without the 
electron-positron contributions, we lack the cancellations that typically take place between 
contributions of the same order. 



If we retain only the self-energy contrit^tions from the two-photon intermediate states, 
the equations for the two-body wave functions become much simpler, and the coupled in- 
tegral equations can be reduced to the one-electron sector. There, they can be solved 
analytically, except for the calculation of certain integrals [13] • We see that the inclusion 
of the self-energy contribution is a significant improvement over the one-photon truncation. 
Thus, we expect that inclusion of the three-photon self-energy will improve the two-photon 
results. 



Figure [5] also includes results obtained for the two-photon truncation when only the one- 
loop chiral constraint is satisfied. Without the full nonperturbative constraint, the results are 



very sensitive to the PV photon mass /ii. This behavior repeats the pattern observed in [12 



for a one-photon truncation without the corresponding one- loop constraint. The resulting 
[i\ dependence is illustrated in Fig. 2 of [l2[. Thus, a successful calculation requires that 
the symmetry of the chiral limit be maintained. 
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FIG. 5. The anomalous moment of the electron in units of the Schwinger term (a/27r) plotted 
versus the PV photon mass, /ii, with the second PV photon mass, //2, set to V^fii and the PV 
electron mass mi equal to 2 • 10^ rUe- The solid squares are the result of the full two-photon 
truncation with the correct, nonperturbative chiral constraint 15|. The open squares come from 
use of a perturbative, one-loop constraint. Results for the one-photon truncation Ij] (solid line) 
and the one-photon truncation with the two-photon self-energy contribution ij] (filled circles) 
are included for comparison. The resolutions used for the two-photon results are = 50 to 150, 
combined with extrapolation to K = oo, and A^_l = 20. 
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